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Detection of a gravitational- wave stochastic background via ground or space-based gravitational- 
wave detectors requires the cross-correlation of the response of two or more independent detectors. 
The cross-correlation involves a frequency-dependent factor — the so-called overlap reduction func- 
tion or Hellings- Downs curve — that depends on the relative geometry of each detector pair: i.e., the 
detector separations and the relative orientation of their antenna patterns (beams). An incorrect 
formulation of this geometrical factor has appeared in the literature, leading to incorrect conclusions 
regarding the sensitivity of proposed detectors to a stochastic gravitational-wave background. To 
rectify these errors and as a reference for future work we provide here a complete, first-principles 
derivation of the overlap reduction function and assess the nature of the errors associated with the 
use of the incorrect expression that has appeared in the literature. We describe the behavior of the 
overlap reduction function in different limiting regimes, and show how the difference between the 
correct and incorrect expressions can be understood physically. 

PACS numbers: 04.80.Nn, 95.55.Ym, 04.30.-w 



I. INTRODUCTION 



The measured response of a single gravitational- wave detector to a stationary stochastic gravitational-wave signal is 
indistinguishable from unidentified instrumental noise. The gravitational-wave contribution to the measured response 
of two or more independent detectors will, however, be correlated between detector pairs in ways that other technical 
noises will not. The relationship between the power in a stochastic gravitational- wave background and the cross- 
correlated response of a detector pair depends on the response of the individual detectors and their relative geometry: 
i.e., their separation and the relative orientation of their respective detector antenna patterns, or beams. In the context 
of ground or space-based laser interferometric detectors [l|, 0, d, 0, H, H, 0] or resonant acoustic gravitational- wave 



detectors @ this geometrical factor is referred to as the overlap reduction function ITol. [Til [l2l. fl3t ; in the context 
of pulsar timing arrays flij or spacecraft doppler tracking [l5j it is called the Hellings-Downs curve 16f. 

Incorrect expressions for the overlap reduction function have appeared in the recent literature [171 [l8| and, with 
them, incorrect conclusions regarding the sensitivity of proposed gravitational wave detectors to stochastic gravita- 
tional waves. These errors have lead to significantly flawed appraisals of the high-frequency sensitivity of the Big 
Bang Observer to a stochastic gravitational wave background, including spurious nulls in the frequency-dependent 
detector response and a reduced estimate of the signal-to-noise ratio as a function of the gravitational- wave power. 
To rectify these errors and as a reference for future work we provide here a complete, first-principles derivation of the 
overlap reduction function and assess the nature and physical interpretation of the errors associated with the use of 
the incorrect expression that has appeared in the literature. 
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II. THE OVERLAP REDUCTION FUNCTION 



The overlap reduction function of a pair of gravitational-wave detectors is the collection of geometric factors, 
associated with the relative position and orientation of the detector pair, that appear in the cross-correlation of 
the detector pair's response. Here we derive an expression for the overlap reduction function by deconstructing the 
cross-correlation, identifying those contributions that depend only on the radiation and those that depend only on 
the detectors, which are then identified as the overlap reduction function. This approach has the virtue of clearly 
illustrating the physical origins of the overlap reduction function and making less likely mistakes of the kind that may 
have led to the existing errors in the literature. 



A. Inter-detector cross-correlation 



Consider two gravitational-wave detectors. A stochastic gravitational-wave "background" will manifest itself in 
a non-vanishing cross-correlation between the measurements unit) made at the two detectors, calculated as a time 
average over the product of the measurements: 

1 f T/2 

C T (At,t) = - dt! m 1 (t + At + t')m 2 {t + t'). (2.1) 

T J -T/2 

The signature of a stochastic gravitational-wave signal hij(t,x) is just the expectation value of Cr(At, t) in the 
presence of hij(t,x). Write the measurement TOj(t) made at detector / as the sum of a noise contribution ni(t) and 
a signal contribution rj(t), corresponding to the detector response to hij(t,x) [23]. Assume that the noise in each 
detector is independent and that there are no non-gravitational-wave effects that might lead to a correlation in the 
measurements made at each detector. Under these assumptions the expectation value of the product ni(t + At)ri2(t) 
vanishes, implying that the expectation value of C<r(Ai, t) is just the expectation value of the product ri(t + At)r2(t): 



C(At) = C T (At,t)=r 1 (t + At)r 2 (t), (2.2) 

where overbar denotes expectation value and we have assumed that the expectation value of Cr(At, t) is independent 
of t. This is equivalent to assuming that the background is stationary, as we will describe in more detail in Sec. Ill CI 
Turn now to the detector response. Gravitational waves are weak. Even the most sensitive detectors respond 
linearly to the local field hij(t,x). Correspondingly we write the detector response as a convolution, in time and 
space, of an impulse response function R l j (t,x) with the field hij(t,x): 

rj(t) = rj(t,xi)= / dr d 3 x hij(t — t,xj — x)R % j {t,x) , (2-3) 
J-oo Jr 3 

where xi is the spatial location of detector / about which its response Rf (t, x) is defined. Causality requires that 
R 1 / (t, x) vanishes outside the future light cone of (0, 0). Exploiting the convolution theorem we can also write 

rj (t) = (2n) 3 ^ df [ d 3 k Taj (/, k)R« (/, fc) e *(^/*-^) (2.4) 
J-oo Jr 3 



where 



Hj (/, k) = ^3 |°° dt d 3 x h tJ (t, Z) e -^n-ls) i (2 . 5a) 



Vf) = 7^3 r dt I ^aitfftaOe-*^*-* 50 (2-5b) 
(2tt) j J_ oc J r3 

are the field Fourier modes and detector transfer function. Note that for real and iZ y 

=M-/,-fc), (2.6a) 

R i i*(f,k)=R i I j (-f,-k). (2.6b) 

With the above representations of the detector response, we can express C(At) in terms of the detector response 
and the field: 

dr / dr' / d 3 x / d 3 x hij(t + At ~ r, 1 1 - x)h a {t -t',x 2 - x')R\ j (T, x),R^ 1 (t',x') (2.7a) 

-oo J-oo JR3 J R? 
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or, equivalently, 



C(At) = (2tt) £ 



df / df 



d A k / <Tk 



•{• 



hij(f,k)h* kl (f>,k>) R iU,k)e 



—ik-. 



2ni(f-f")t p 2nifAt 



(2.7b) 



Note particularly how the detector location and the transfer function appear together in the combination 
Ftj (/, k)e~ lk ' Xl . The form of this combination will be critical when we come to understand the physical charac- 
ter of the errors made in earlier calculations of the overlap reduction function. 



B. Plane- wave representation of stochastic signal 

Focus attention on gravitational wave fields hij(t, x). These are conveniently represented as a superposition of plane 
waves 



oo JS 2 



ha(t, x) = df d 2 n~ k e 2 ™M-~ k ^H A (f, k)e?j(k) , 



(2.8) 



where k is the unit vector direction of wave propagation, and e A j(k) are the two orthogonal polarization tensors, 



26 AA =ef j (k)ef j (k) , 

A 

ji 



et i (k)=ef i (k) 



(2.9a) 
(2.9b) 



Note that H A (f, k) = HA(—f, k) as a consequence of the reality of hij. The plane-wave field amplitudes Wa(/, k) are 
related to the field's Fourier modes hij (/, k) by 



H A (f,k)eUk) 



(2vr/) 2 H|+ ) (/,fc), />0 
(2vr/) 2 H l ( 7 ) (/,-fc), /<0 



(2.10) 



where 



hij(f,k) = H\+\f,k)S(\k\-2nf) + H\-\f,k)5(\k\+2nf). 



(2.11) 



Here we have introduced separate amplitudes for the positive and negative frequency solutions to the dispersion 

relations \k\ 2 = (27r/) 2 for a plane wave. 

Using expansion (|2.8[) . we can write the detector response rj(t) as 



oo J S 



ri(t) = I df I d^HAif^mfif,^ 2 "^-'^ 



where 



(2.12) 



Kf{f 1 k) = (2vfetAk)R l I \f,2vfk) 



(2.13) 



The expectation value C(At) can also be written as 



c(At) = f df f df f d 2 n k f d 2 n~ k , {, 

J-oo J-oc JS 2 JS 2 1 



2Ki(f-f')t27TifAt 



n A (f,k)n A ,(f',k') [izf(f,k)e- 2 ^' k -^] [n A \f,k')e~ 2 ^f' y - S2 \* } 



(2.14) 



Again make note of how the detector locations x\, X2 are associated with the respective transfer functions TZ A and 



K A . 
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C. Stationarity, isotropy, and polarization correlations 

The statistical properties of the stochastic signal are encoded in the expectation values of products of the gravita- 
tional field 



hij(t,x), hij{t,x)h k i{t',x') , hij(t, x)h k i(t',x')h mn (t",x") , • • ■ (2.15) 

Without loss of generality, we will assume that any non-zero mean has been absorbed in the background spacetime, 
so that hij(t, x) = 0. Furthermore, for Gaussian-distributed fields, knowledge of the quadratic correlations will suffice 
as all higher-order moments can be constructed from these. 

In our problem we expect that the gravitational wave background is effectively stationary: i.e., that (t, x)hu (f , x') 
depends on t and t' only through their difference t — t'. In terms of the plane wave components TL A (f, k), this condition 
becomes 



H A (f, k)H* A , (/', k') = H AA , (/, k, k')6(f - /') . (2.16) 

Thus, the different frequency components of a stationary stochastic background are statistically independent, but they 
can contribute differently to the cross-correlated power through the /-dependence in H^/. Note that the 5(f — /') 
factor in (|2.16|) eliminates the i-dependence in C(At), cf., Eq. 12.141 

If the background is isotropic — i.e., the gravitational- wave specific intensity is independent of the direction of 
propagation k — then the most general form of the quadratic expectation value of the plane wave components 

n A (f,k) is 



H A (f, k)TL* A , (/', k') = Haa'U, /', k ■ k') , (2.17a) 

where H A A' depends on k and k! only through the angle between them. If we further assume that the components 
corresponding to different propagation directions are statistically independent, then 



n A (f,k)n%(f,k') = H^(/,/')<5 2 (fc,fc'), (2.17b) 

where 5 2 (k, k') = S(cos9 — cos — cfi') is the covariant Dirac delta function on the two-sphere. This latter, more 
restrictive, condition is the definition of isotropy for gravitational-wave stochastic backgrounds typically assumed in 
the literature, e.g., (TSUI H0J] • 

Finally, if the background is unpolarised, by which we will mean that the different polarisation components are 
statistically independent and contribute equally to the cross-correlated power, then 



H A (f, k)n* A ,(f, k') = H(/, /', k, k!)8 AA , . (2.18) 

Putting all these conditions together, we have that an unpolarised, stationary, isotropic stochastic gravitational-wave 
background satisfies 



n A (f, k)H* A ,(f, k') = H(f)S(f - /> 2 (fc, k')5 AAI . (2.19) 

Here H(/) is a real- valued function proportional to the gravitational- wave energy density. This is, in t urn, directly 
related to fi gw (/), the ratio of the gravitational- wave energy density to the cosmological closure density |13| : 

"W-iSW- (2 - 20) 

where Hq is the Hubble expansion rate at the present epoch. 



D. The overlap reduction function 

Combining the results of the previous subsections we find that we can express the expectation value C(Ai) of the 
inter-detector cross-correlation in the presence of an unpolarised, stationary, isotropic gravitational-wave background 

C(At)= I d/e 2 " /At H(/)r 12 (/), (2.21a) 
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where 

r 12 (/) = / d 2 n k K?(f, k)1Z£*(f, fc) e - 2 «^-(^-^) , (2.21b) 



s 

oo 



Kf(f,k)=ef j (k) [ dt f dfxRWt.xy- 2 * 1 ^ 1 ^ . (2.21c) 

J-oo J R 3 



The quantity Ti2(/) is the overlap reduction function. It is often convenient to define a normalized overlap reduc- 
tion function 712 (/) oc T^/) with 712(0) = 1 for two coincident and coaligned identical detectors. For identical 
interferometers with opening angle (3 this leads to the normalized form 
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8tt sin" p~ 



JM = n _^ 2 j Mf). (2.22) 



Some symmetry properties of 712 (/) follow from immediately from its definition: in particular, 

7i* 2 (/)=72i(/), (2.23a) 
7i* 2 (/)=7i2(-/). (2.23b) 



III. DISCUSSION 



Summarizing the results of the previous section, the overlap reduction function normalized for interferometric 
detectors is 

112(f) = — ^r-r I d 2 n k n?(f, k)n$*(f, %) e -™*<*-** , (a.ia) 

87r sin p Js 2 

where 

/OC l> 
dt d 3 a;^(t,f)e- 2 " /( *-^ s) . (3.1b) 
-00 Jr 3 

Adjusting where necessary for differences in the plane- wave expansion and detector numbering conventions, references 
[l7l [l8| give the overlap reduction function as 

7i 2 (/) = 7T^-r I d 2 n k nt(f,k)Tlt(f,k)e 2 ^' k < s ^K (3.2) 
»7r sin p Js 2 

Comparing this expression with our Eq. I3.1al we see that they differ by the substitution of e 27T tfk-{x 1 -x 2 ) f or 
e -2TTifk-(x 1 -x 2 ) _ xhis difference and when it is significant can be understood physically; doing so provides the occasion 
for a deeper discussion of the overlap reduction function. 



A. Detector locations 



At the end of Sec. Ill Bl we observed that the response of detector /, located at xi, to a field of plane gravitational 
waves is given by (cf. Eq. 12. 12ft 

ri (t) = r df [ d 2 n k n A (f, k)nf(f, k) e 2 -/c-^) . (3.3) 

Js 2 

The detector location appears here in the form e -21 "^' 1 '. Referring to Eqs. I3.1al and I3T21 for 712 and 7^ it is clear 
the substitution of 7 12 for 712 is equivalent to simply exchanging the locations of detectors 1 and 2 keeping the rest 
of the configuration of the detectors fixed: i.e., 7(2 is the overlap reduction function for the detector configuration 
consisting of detector 1 at location xi and detector 2 at location x\. With this understanding we now ask when that 
exchange is significant and when it is not. 
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B. Radiation wavelength and detector separation 

An intuitive understanding of 712 (/) recognizes that its behavior in different frequency regimes is governed by 
several independent dimensionless parameters that can be created from the radiation wavelength, the separation 
between the detectors, and several intrinsic properties of the detectors as they are represented in the detector impulse 
response functions. 

Referring to Eq. 13. lal we note that when f\xi — x 2 \ <C 1 the exponential term may be replaced by unity. As this is 
the only place where the detector separation appears, in this limit the detector separation plays no role in determining 
the value or behavior of "/12(f)- Defining 

5 = f\xi-x 2 \ (3.4) 

we refer to 8 <C 1 as the small separation limit. In the small separation limit, then, the difference between 7 12 and 
7^2 is negligible under all circumstances. 

Now consider the case 6 > 1. As we have observed, the difference between 712 and 7(2 is the difference between 
locating detector 1 at x\ or x 2 , an d detector 2 at x 2 or x±. When the two detectors are identical in all other aspects, 
so that IZi — 7I2, this exchange leaves the physical configuration unchanged and, again, there will be no difference 
between 712 and 7 12 . 

To understand the case 8 > 1 when the two detectors are not identical we must consider the detector impulse 
response functions as they appear in Eqs. 13. ll 

C. Radiation wavelength and detector impulse response 

The impulse response of a detector has finite support: i.e., R l /(t,x) ~ for sufficiently large t > or \x\. A detector 
doesn't sample the field beyond its physical extent, so the support in x will be on order the detector's size I. Referring 
to Eq. I3.1bl it is apparent that lZf(f,k) depends on k only through e^(fc) when f£ <C 1: 

^(/,fc)^(2^) 3 e^(fc)^'(/,0). (3.5) 

Introducing the parameters 

6! = fll (3.6) 

for the two detectors 7=1,2 we refer to ej <C 1 as the small antenna limit for detector I. 

Return now to the difference between 712 and 7 12 when 8 > 1 and the detectors are not identical. When both e\ 
and 62 are <C 1 and noting that efj(-k) = ±e^(fc) we find 

7l2 (/) ~ — (2tt) 6 (/, 0)i?^ * (/, 0) / d 2 ^e^(fc)e^(fc)e- 2 ^-^-^' 
87rsm p Js 2 

= ^T^(2^) 6 i?r'(/,0)i? 2 fci *(/,0) / d 2 0_ fe e4(-fc)e^(-fc) e 2 ^-^-^) 
87rsm p J S 2 

= ^T^(2^) 6 i?r'(/,0)i? 2 fci *(/,0) / d 2 ^e^.(fc)e^(fc)e 2 ^-(---)^7i 2 (/); (3-7) 
87rsm p Js 2 

i.e., in the small antenna limit there is no distinction between 712 and 7 12 . 

D. Large separations and large detectors 

Finally, consider the case 8 > 1 and, without loss of generality, e% > 1. In this case, the variation of the field across 
the spatial extent of detector 1 is important to the detector response and, in turn to 712. Exchanging the detector 
locations changes the relationship between the spatial extent and orientation of detector 1 relative to the location of 
detector 2. Correspondingly, in this limit the distinction between 712 and 7 12 is significant. 

An example of the case 8 > 1, e\, e% > 1 is the anticipated sensitivity to a stochastic gravitational- wave background 
of the Big Bang Observer (BBO), a space-based follow-on to LISA that has been the subject of recent study 0, [l?], 
EH, H3j El ■ The principal results of this study reported in the literature [l7], [H| make use of the incorrect form of the 
overlap reduction function, thus mis- estimating this proposed detector's sensitivity to a stochastic gravitational- wave 




Detector 1 



Detector 2 




FIG. 1: The arrangement of the two BBO detectors described [13] and the arrangement actually analyzed when y[ 2 is substituted 
for 712. Owing to its larger spatial extent, the system actually studied is much more sensitive to frequency than the system 
whose study was intended. 



signal in the higher frequency regime. Figure [T] illustrates the the physical effect of using y' l2 m place of 712 when 
analyzing the cross-correlation of the two BBO detectors as described in [13] ■ On the left are the two interferometric 
detectors as they are actually arranged in space; on the right are the effective location and orientation of the detectors 
when 7^2 is used in place of 712: i.e., when detector 1 is translated to x-i and detector 2 is translated to X\. Under this 
transformation the spatial extent of the two detector pairs is much greater than is actually the case; correspondingly, 
we expect j' 12 to be a much more sensitive function of frequency than 712. 

To calculate 7(2 for comparison with 712, we need an explicit expression for the transfer function of the detectors. 
This is derived in (l7l Eqs. 5, 7, 11]. In our notation, 



TZ A (f, k) = eM)I ( t*VT(fi • k, /) - «VT(fi • k, /) 



(3.8a) 



where 



T(u ■ k, J) 



smc 



/ 



(1 — u ■ k) I exp 



smc 



(1 + u • k) 1 exp 



/ 



(1 + u-k) 



(3.8b) 

Here u and v are unit vectors pointing in the direction of the detector arms, /* = c/(2ttL) is the transfer frequency 
of the detectors (arm length L = 5 x 10 9 m), and sinc(x) = sin(a;)/x. Figure [2] compares 7 12 (cf. [U Fig. 5]) and 712 
(Eq. 13. lj) for this configuration. As expected 7^2 decreases much more rapidly with frequency and has its nulls more 
closely spaced than those of 712. Figure [3] shows the fractional error |1 — 7i2(/)/7ia(/)| as a function of frequency. 
The large amplitude spikes in the error occur where nulls of "/' 12 do not coincide with nulls of 712 . 
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10" 2 1Cf 1 „,„ 10° 10 1 



FIG. 2: Plots of the overlap reduction functions 712 (/) (solid line) and 712 (/) (dot-dashed line) for the BBO configurations in 

Fig.ru 
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It is apparent from Fig. [2] that the use of an incorrect overlap reduction function has, in this case, led to an 
underestimate of BBO's sensitivity to a stochastic gravitational wave background. Estimating the background and 
detector noise power spectral density as white, the magnitude of the error is just the ratio of integrated squared 
magnitudes |7i2(/)| 2 and |7i 2 (/)| 2 , 



(3.9) 
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i.e., the BBO estimates of [13, EH underestimate the sensitivity of BBO by nearly 30% across the entire band, 
and substantially larger if interest is focused on the higher frequencies. Achieving BBO's goals of detecting the 
stochastic gravitational-wave relics of the inflationary epoch depend on the accurate identification and subtraction of 
contributions owing to compact binary systems. Underestimating BBO's response to a gravitational- wave background 
leads to an overestimate of the accuracy required in this identification and subtraction [19( . Recognizing and correcting 
the underestimate in BBO sensitivity thus relaxes the analysis problem associated with the identification of these 
foreground sources. 

IV. CONCLUSIONS 

Detection of a gravitational- wave stochastic background relies on the cross-correlated response of one or more pairs 
of gravitational- wave detectors. The separation and relative orientation of the two detectors plays a crucial role in 
determining the frequency dependent sensitivity of each detector pair to the stochastic background. Recent studies 
[13, [HI of the sensitivity of the Big Bang Observer and related future generation gravitational- wave detectors have used 
an incorrect expression for this geometrical factor. The errors committed may be physically interpreted as an exchange 
in space of the two detectors, leaving their absolute orientations fixed. In the case of the Big Bang Observer, this error 
leads to an approximately 30% underestimate in its sensitivity to relic gravitational waves associated with, e.g., the 
inflationary epoch. Since achieving BBOs goals of detecting this background requires the accurate identification and 
subtraction of gravitational- wave foreground contributions from compact binary systems [l9T ] , this underestimate has 
lead to commensurate overestimate of the difficulty of this analysis problem. Recognizing and correcting this error 
thus improves, in two ways, the prospects for the BBO missions main goal as a Big-Bang Observer. 
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